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Based on the full band electronic structure calculations, first we consider the effect of n-type 
doping on the optical absorption and the refractive index in wurtzite InN and GaN. We identify 
quite different dielectric response in either case; while InN shows a significant shift in the absorption 
edge due to n-type doping, this is masked for GaN due to efficient cancellation of the Burstein-Moss 
effect by the band gap renormalization. Moreover, for high doping levels the intraband absorption 
becomes significant in InN. For energies below 1 eV, the corresponding shifts in the real parts of 
the dielectric function for InN and GaN are in opposite directions. Furthermore, we observe that 
the free-carrier plasma contribution to refractive index change becomes more important than both 
band filling and the band gap renormalization for electron densities above 10 19 cm -3 in GaN, and 
10 20 cm -3 in InN. As a result of the two different characteristics mentioned above, the overall change 
in the refractive index due to n-type doping is much higher in InN compared to GaN, which in the 
former exceeds 4% for a doping of 10 19 cm -3 at 1.55 /im wavelength. Finally, we consider intrinsic 
InN under strong photoexcitation which introduces equal density of electron and holes thermalized 
to their respective band edges. The change in the refractive index at 1.55 um is observed to be 
similar to the n-doped case up to a carrier density of 10 20 cm -3 . However, in the photoexcited 
case this is now accompanied by a strong absorption in this wavelength region due to Y% —¥ 
intravalence band transition. Our findings suggest that the alloy composition of In^Gai^N can 
be optimized in the indium-rich region so as to benefit from high carrier-induced refractive index 
change while operating in the transparency region to minimize the losses. These can have direct 
implications for InN- containing optical phase modulators and lasers. 



PACS numbers: 78.20.Ci, 78.20.Bh, 78.40.Fy 



I. INTRODUCTION 



The research efforts on InN and In-rich InGaN have 
been intensifying worldwide. As well as their traditional 
applications on solid-state lightning and lasers, new pos- 
sibilities begin to flourish such as photovoltaics and chem- 
ical sensing^ From the electronic structure point of view, 
InN has a number of unique properties^ Its band gap for 
the wurtzite phase is commensurate with the 1.55 /Ltm 
wavelength for fiber optics when it forms an alloy with 
a small amount of GaN. Furthermore, the small conduc- 
tion band effective mass enables interesting band-filling 
effects. Even though the latter giving rise to Burstein- 
Moss effect 3 -^ was the center of focus in the recent re- 
assessment of the high quality InN samples^— its device 
implications have not been given the full attention it de- 
served. Among the viable applications are the carrier- 
induced optical phase modulators^ and tunable Bragg 
reflector and filters i 10 ' 11 These devices are expected to 
play essential role, for instance, for the signal processing 
directly in the optical domain in the future high bit rate 
optical communication networks. 

It was known from the early days of edge-emitting 
semiconductor lasers that carrier-induced refractive 
change affects the optical beam quality^ such as 
mode guiding along the junction plane, self-focusing, 
filament formation, frequency chirping under direct 



modulation! 13 ' 14 The refractive index tunability by car- 
rier injection was studied in the previous decade for 
the Indium-Group-V semiconductor compounds and al- 
loys other than InNji£~— In the case of InGaN laser 
structures, there have been recent theoretical^ and 
experimenta l 19 ' 20 studies on the carrier-induced refrac- 
tive index change and the linewidth enhancement factor. 
However, these studies considered only Ga-rich limit of 
InGaN alloys. 

The purpose of this work is to offer a comparative theo- 
retical account of carrier-induced refractive index change 
in both InN and GaN, so as to form a basis for the In- 
rich InGaN alloys. Primary attention is given to n-type 
doping which is the prevalent type among the routinely 
grown samples.— For the case of InN, the photoexcited 
bipolar carrier-induced dielectric effects are also exam- 
ined. The crux of our study is based on the Bennett- 
Soref-Del Alamo approach-^ which includes band-filling, 
band gap renormalization (BGR) and plasma contribu- 
tions. However, rather than a model band structure, we 
employ a rigorous full band technique to represent the 
band filling; preliminary account of this work was given 
in Ref. 21. In Sec. II we present the theoretical details, 
this is followed by our results in Sec. III. A self-critique of 
our model and further discussions are provided in Sec. IV, 
ending with our conclusions in Sec. V. The Appendix con- 
tains our empirical pseudopotential parameters used for 
wurtzite InN. 



II. THEORY 

First, we would like to describe the electronic struc- 
ture on which our computations are based. In the case of 
GaN, we use the local empirical pseudopotential method 
that we fitted to the existing experimental data, putting 
special emphasis on the conduction band behavior. 22 
For InN, we use the nonlocal empirical pseudopotential 
method, treating the angular-momentum channels p and 
d for In, and p for N as nonlocal. 23 Further details are 
given in the Appendix section. The corresponding band 
gap for InN comes out as 0.85 eV which agrees very well 
with the recent first principles result o 24 ' 25 but it is higher 
than the 0.64 cV value established experimentally. 1 Fig- 
ureQ]shows the band structures of wurtzite InN and GaN. 

As the Burstein-Moss effect 3 -" 4 " primarily shifts the ab- 
sorption edge due to band filling, the ultimate quantity 
to be computed is the imaginary part of the optical (i.e., 
the long wavelength limit) dielectric tensor which is given 
in Gaussian units by 2 ^ 

Im { e ab H} = - E / dkr a vc (k)r b cv (k) 
n v,c Jbz 

xS(E c (k) -E v {k) -Hu) , (1) 

where a, b are the Cartesian indices. r^ c ik) = 
Py C {k) I (imoU} vc (k)), where p^ c (k) is the momentum ma- 
trix element, too is the free electron mass, uj vc (k) = 
uj v (k) — uj c (k), where Huj n (k) = E n (k) is the energy of 
the band n, at the wave vector k. The specific band la- 
bel v (c) represents filled (unfilled) bands, where for the 
clean cold semiconductors, the band occupation factors 
are taken as either one or zero. The volume integration 
in Eq. (JT|) is over the first Brillouin zone (BZ) which can 
be reduced to the irreducible BZ using the symmetry re- 
lations of the BZ (see, Fig. 
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FIG. 1. The empirical pseudopotential band structures for 
InN (solid) and GaN (dashed) used in this work. 



FIG. 2. Schematic (not to scale) of the Brillouin zone of the 
wurtzite structure (top), its irreducible Brillouin zone (mid- 
dle), and its tessellation (bottom); each volume is further di- 
vided into tetrahedra (not shown), a and c are the lattice 
constants in the basal plane and perpendicular directions, re- 
spectively. 



The Dirac delta term in Eq. ([l) automatically re- 
duces the volume integration to a surface formed by 
the k points which allow direct (vertical) transitions 
from a filled state to an unfilled state with an en- 
ergy difference corresponding to the chosen photon en- 
ergy. Such surface integrations routinely appear in the 
density of states, effective mass and response function 
calculations; 27 They can be efficiently calculated numer- 
ically using the Lehmann-Taut method, 28 This technique 
has also been shown to work for the second-order non- 
linear optical response function 2 ^ and also for the full 
band phonon scattering rate calculation^ In this work, 
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for the required accuracy, we divide the irreducible BZ 
into a mesh of 40 x 40 x 40 along the basal plane and 
the c-axis directions (see, Fig. [2J. In Fig. we show 
the imaginary part of the optical dielectric function for 
the intrinsic (i.e., undoped) InN and GaN. The ordinary 
tensor component which corresponds to the electric field 
parallel to basal plane of the hexagonal crystal (i.e., in- 
plane) polarization is shown. For GaN, we observe a 
discrepancy between the theoretical and experimental 
values close to the band edge. On the other hand, our 
theoretical result for InN matches well with the rigorous 
first-principles studies; 32 ' 33 and also in reasonable agree- 
ment with the experiment £Q 
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FIG. 3. Imaginary part of the in-plane optical dielectric func- 
tion for the intrinsic InN (blue) and GaN (red), comparing 
our theoretical (solid) results with the experimental (dashed) 
data of Refs. [30 and[3ll. 



in the case of InN. Therefore, in these expressions we 
use an energy-dependent density of states effective mass 
value corresponding to the level of band- filling. Figure |4] 
shows for the conduction band these values used in this 
work for InN and GaN. The curve for InN is taken from 
the first-principles results^ For GaN, it has been com- 
puted using the same Lehmann-Taut method, but with 
a finer grid than used above, such as 60 x 60 x 60 or 
more. In the case of photoexcitation which we consider 
for InN, actually we have a two-component quantum sys- 
tem (i.e., electron and hole gases). In the interests of sim- 
plicity, we adapt Eq. © by replacing the electron effec- 
tive mass with the reduced electron- hole effective mass, 
Heh given by — m" 1 +m~ h 1 where m e and m/, are the 
energy-dependent electron and hole effective masses, re- 
spectively. The latter is almost constant for the densities 
dealt in this work and therefore is taken to be 0.84mrj for 
InN. 



0.40 -| 



0.35- 




0.05 I 

10 17 10 18 10" 10 20 10 21 

Electron Concentration (cm 3 ) 



Another physical mechanism that accompanies band- 
filling is the BGR due to many-body interactions among 
the carriers* 3 ^ In highly n-doped InN samples it is 
observed that this partially cancels the Burstein-Moss 
effect.' The electron-electron and electron-ion contribu- 
tions within the random phase approximation are given 
in Gaussian units byZiM 



AE e - e = 



e 2 kxF 



2c 



4 f k F 

1 arctan 

7T V ktf 



Ane 2 N 



(2) 
(3) 



where N is the electron (ion) density and e s is the static 
permittivity for which we use the experimental value 
of 9.5 (15.3) for GaN (InN)^^ k F = (3tt 2 N)^ 3 and 

kxF = 2yfe^7(™fl) are ' respectively, the Fermi and 
Thomas- Fermi wave numbers, a* B = 0.529 e s mo/m* is 
the effective Bohr radius in A, m* is the conduction band 
effective mass. The latter quantity rapidly deviates with 
increase of energy from its band edge value, especially 



FIG. 4. The density-of-states effective mass variation as a 
function of band-filling for the conduction band electrons in 
InN and GaN. The latter is computed in this work, whereas 
the former is taken from Ref. Y2A. 



Having determined the imaginary dielectric function 
of a semiconductor with partially filled band(s), sub- 
ject to BGR, the next aim is to compute their effect on 
the refractive index. Their interrelation is established 
by the causality principle which has to be satisfied by 
any physical response function^ In mathematical terms, 
this is governed by the Kramers-Kronig relation; for the 
real part, the corresponding expression, suppressing the 
Cartesian indices is given by 



Re{e(uj)} = l + -V [ 
7T Jn 



(4) 



Computationally, it requires the knowledge of the imag- 
inary part of the dielectric function for all energies; we 
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have verified the convergence of our results by including 
60 bands and accounting for the higher energy behavior 
(above 30 eV) analytically with a 1/uj 2 fall off. 

Based on the mechanisms described above, carrier in- 
jection essentially leads to a change in the refractive in- 
dex. Another significant contribution to the refractive 
index change is caused by the plasma absorption of the 
free carriers: electrons in the conduction band for the 
n-doped case or the electrons and holes for the photoex- 
cited case. Its resultant contribution can be expressed in 
convenient units as 15 



6.9 x 10" 



plasma 



N[c 



n m* E 2 [eV] 



(5) 



where no is the refractive index of the intrinsic semicon- 
ductor, E[eV] is the energy in eV, A[cm~ 3 ] is the free- 
carrier density in cm -3 and m* is the effective mass in 
units of mo; in the photoexcited case, for the latter two 
we use the electron-hole density and the effective mass 
(peh), respectively. 



III. RESULTS 
A. n-doped case in InN and GaN 

We first start with the n-doped case which is quite eas- 
ily achieved (sometimes even unintentionally) in InN and 
GaN and, we explore the filling of the conduction band by 
electrons fully ionized from their dopants. The electron 
concentrations up to about 7xl0 20 cm -3 (9xl0 20 cm~ 3 ) 
are considered marking the threshold beyond which the 
next higher conduction band starts to be filled for InN 
(GaN). Figure [5] shows the shift in the absorption edge 
due to Pauli blocking of the filled portions of the lowest 
conduction band. A marked contrast revealed with this 
plot is that the absorption edge shift in InN is much 
more significant compared to GaN. The disparity be- 
tween the two materials' conduction band edge density 
of state behavior plays a role here. While in InN the 
band-filling effect is quite strong, in the case of GaN this 
becomes smaller which is moreover effectively canceled 
by the BGR. 

Yoshikawa et al. have used temperature-dependent 
photoluminescence spectroscopy to separate BGR and 
band-filling effects in Si doped GaNj2& Figure[S]iHustrates 
the comparison of our n-doped GaN results with this 
work for the reduced band gap, i.e., due to BGR, as well 
as for the optical gap which corresponds to the absorption 
edge including the BGR and Burstein-Moss effects. The 
dashed line is the fitted expression, — 4.72 x lO^n 1 / 3 eV 
offered by Yoshikawa et al. for the former. Very re- 
cently Schenk et al. relying on band edge luminescence 
peak and the line shape information have extracted an 
expression given by —2.6 x 10~ s n}/ 3 eV for the BGR 
in n-doped GaN* 3 -^ While our result closely agrees with 
the Yoshikawa data, Schenk's expression is in disagree- 
ment with these. In Fig. [51 also a much more elaborate 



theoretical analysis for the BGR is shown which how- 
ever deviates from the rest of these results^ We believe 
that further studies are required to reach to a consensus 
over the BGR expression for GaN. The situation is more 
harmonious for InN. The comparison of our computed 
absorption shift for InN with the measurements and cal- 
culation from Wu et al.^ together with other data men- 
tioned in that worker— is provided in Fig. [7] it shows 
the overall agreement of our results with the literature. 
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FIG. 5. The shift in the absorption edge by n-type doping for 
GaN and InN due to conduction band filling without (dashed) 
and with (solid) BGR accounted. The inset shows the density 
of states (DOS) for InN together with the band-filling levels 
on the valence and conduction bands in the absence of BGR. 



The conduction band filling and the associated BGR 
directly affect the imaginary part of the dielectric func- 
tions through the shift in the absorption edge deeper 
into the conduction band. This is illustrated for the or- 
dinary components (i.e., the electric field lying on the 
crystal basal plane) for different electron concentrations 
in Figs. [5] and for InN and GaN, respectively. It can 
be noticed for both cases, especially for larger dopings, 
that there is a nonzero absorption below the absorption 
edge, even down to zero energy. This is due to intraband 
absorption within the conduction band (see the upper 
inset in Fig. [5]). The same behavior also appears in the 
measured InN samplesi* 

The real parts of the dielectric functions are obtained 
from the imaginary parts using the Kramers-Kronig rela- 
tion in Eq. The results are shown in Figs. [TU1 and [TT1 
for InN and GaN, respectively. First of all, the real part 
of dielectric function of InN is more sensitive to n-type 
doping. Furthermore, for low energies, the shifts are in 
opposite directions for InN and GaN, so that as the elec- 
tron concentration increases, the permittivity decreases 
for InN, whereas it increases for GaN. Note that these 
are only the electronic contributions to polarization as 
Eq. ([I]) does not include the ionic degrees of freedom. 
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FIG. 6. The comparison for the n-doped GaN, the reduced 
gap and the optical gap values (see text) with Yoshikawa el 
aL,— Schenk el aZ.— and Persson el aZ.— The band gap of 
intrinsic GaN is taken as 3.5 eV. 



FIG. 8. The calculated imaginary part of the in-plane dielec- 
tric function (DF) of InN for different n-type doping densities. 
Lower inset shows the same graph for a wider energy scale. 
Upper inset indicates the responsible intraconduction band 
transitions for the low energy part of the spectra. 
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FIG. 7. The comparison of measured and calculated shifts in 
the absorption edge in InN with the calculation and measure- 
ments from Wu, 8 together with other data mentioned in that 
work: Haddad,— Inushima,— Tyagai.— 



FIG. 9. The calculated imaginary part of the in-plane dielec- 
tric function (DF) of GaN for different n-type doping densi- 
ties. Inset shows the same graph for a wider energy scale. 



Nevertheless the ions do not directly play a role in the 
change in the refractive index which is focus of this work. 
As mentioned in the Theory section, one should also add 
the free carrier plasma contribution to the refractive in- 
dex [see, Eq. dSJ]. This is separately shown in Fig. [12] 
which always has a negative contribution, and | An p i aS ma| 
decreases as 1/E 2 with energy E. 

The overall change in the refractive index, including 



all of the effects mentioned so far is shown in Fig. [13] 
Because of their technological importance, it is evalu- 
ated at the two fiber-optic communication wavelengths 
of 1.55 and 1.3 /an which correspond to energies of 0.8 
and 0.954 eV, respectively; the former is marked with 
dashed vertical arrows in Figs. [TU] and [TT] It can be 
observed that the refractive index change due to n-type 
doping is much higher in InN compared to GaN for both 
wavelengths. There are several reasons behind this out- 
come. First of all there is a cancellation that occurs for 
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FIG. 10. The calculated real part of the in-plane electronic 
dielectric function (DF) of InN for different n-type doping 
densities. The vertical red arrow marks the 0.8 eV (1.55 um) 
value. Inset shows the same graph for a wider energy scale. 



ally dominates beyond a density of 10 cm' J in GaN 
and 10 20 cm~ 3 in InN. 

Note that at the 1.3 /im wavelength the intrinsic as 
well as doped InN display substantial interband absorp- 
tion (see, Fig. [8]), hence, in Fig.[l3]we used the following 
generalization of the refractive index in the presence of 
loss££ 



v/Re{e} 2 + Im{e} 2 + Re{e} 



(6) 



Because of this loss, 1.3 urn case shows an overall reduc- 
tion in the refractive index change compared to 1.55 um 
wavelength. As another curious point, the n-doped InN 
curve in the bottom panel of Fig. [T3] starts from a pos- 
itive value which requires an explanation. Due to the 
absorption in the 1.3 /jm wavelength, the refractive in- 
dex also peaks around the 0.95 eV range, particularly for 
the 10 18 cm~ 3 density (cf. Fig. [TU|) . Hence, the n-doped 
InN curve in Fig. [13] becomes positive for this density as 
its refractive index exceeds that of the intrinsic sample. 
Even though the band-filling and plasma contributions 
try to reduce it, they are not as effective at this density. 
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FIG. 11. The calculated real part of the in-plane electronic 
dielectric function (DF) of GaN for different n-type doping 
densities. The vertical red arrow marks the 0.8 eV (1.55 /im) 
value. Inset shows the same graph for a wider energy scale. 
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FIG. 12. The comparison of the free-carrier-induced refrac- 
tive index change, i.e., plasma contribution, [see, Eq. (0] for 
two different n-type doping densities for InN and GaN. The 
photoexcited N = P = 10 19 cm -3 bipolar case in InN is also 
included. 



GaN; as shown in Fig. [5j BGR very effectively cancels 
the Burstein-Moss shift up to a density of 10 19 cm -3 . 
The lack of such a cancellation in InN is because of the 
very strong band-filling effect as a result of small density 
of states of InN close to conduction band edge, see the 
inset in Fig.[5j Furthermore, contrary to the case in InN, 
the refractive index change in GaN due to band filling 
is positive up to the ultraviolet. This is canceled by the 
negative plasma contribution [cf., Eq. (|B"])]. which actu- 



B. Photoexcited case in InN 

We would like to compare the pronounced refractive 
index change in n-doped InN with the case where both 
electrons and holes are present. Since introducing holes 
with p-type doping is still a major technical challenge 
in InN,— we consider photodoping 44 where equal density 
of electrons and holes are generated by a strong pho- 
toexcitation. We further assume that these photocarri- 
ers quickly thermalize to their respective band edges and 
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applications. 



GaN (n-doped) 




10" 10 a 
Free Carrier Concentration (cm 3 ) 



10' 



FIG. 13. The percentage refractive index change with respect 
to intrinsic (i.e., no free carrier) case for InN and GaN. All 
effects considered in this work are included, i.e., band-filling, 
BGR, and plasma contributions. Top panel shows the 1.55 /im 
wavelength (0.8 eV) and the bottom panel corresponds to the 
1.3 (jm wavelength (0.954 eV). See text, for the explanation 
of why n-doped InN curve in the bottom panel starts from a 
positive value. The lines are to guide the eyes. 
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FIG. 14. The calculated imaginary part of the in-plane di- 
electric function of InN for different electron-hole popula- 
tions forming degenerate quasiequilibrium distributions. In- 
sets show the intravalence band transitions that give rise to 
the new absorption channels around eV (left) and 0.94 eV 
(right). 



attain a degenerate quasiequilibrium distribution. Since 
the electron-hole pairs will eventually recombine within 
a time scale of 300-400 ps^ the quasicontinuous optical 
excitation should replenish this in order to sustain the as- 
sumed steady-state carrier distribution. The associated 
quasi-Fermi levels for different electron-hole densities are 
marked in the inset of Fig. [5] which clearly indicates that 
the contribution of the holes to the Burstein-Moss effect 
is negligible. In Fig. [14] we show the imaginary dielectric 
function at several electron-hole densities. Compared to 
the n-doped case in Fig. [5] we observe two new features 
around and 0.94 eV which are caused by the intrava- 
lence band transitions: the former is due to very low 
energy transitions at the valence band maximum (left 
inset in Fig. [14]), whereas the latter is the result of the 
Fg -> Tg transition (right inset in Fig. [13]). It is the latter 
that makes a significant impact around 1.55 /im (0.8 eV) 
and even more so at the 1.3 urn wavelength (0.954 eV). 
Therefore, we would like to assert that our band structure 
agrees quite well with the other first-principle studie a 25 i 32 
which report the value 0.9 eV with the quasiparticle cor- 
rection for the Tg — > Tg transition. As observed from 
Fig. [M] these new intravalence band channels refill the 
part of the spectra cleared by the Pauli blocking, undo- 
ing the Burstein-Moss effect. Also we can observe from 
Fig. [T3] that below a carrier density of 10 20 cm~ 3 , the 
change in the refractive index at 1.55 fim. is seen to be 
very close to the n-doped case. However, it should be re- 
called that this is now accompanied by a strong absorp- 
tion (cf., Fig. [H]) which can hamper its possible device 



IV. A SELF-CRITIQUE OF THE 
THEORETICAL MODEL AND DISCUSSIONS 

In this section we would like to discuss several aspects 
of the employed theoretical model both for the assertion 
of its general validity as well as for pointing out its pos- 
sible improvements. To start with, all our results follow 
from the dielectric function expression of Eq. ([1]) which 
is obtained within the independent particle (i.e., random 
phase) approximation^ and as such, it neglects the ex- 
citonic effects. In the case of intrinsic InN dielectric func- 
tion, it has been shown by Furthmiiller et al. that, apart 
from small shifts in the energy positions, these excitonic 
effects are responsible for the enhancement of some of the 
peaks in the spectra, predominantly close to the band 
edge, while the general shape remaining as unaffected. 32 
Furthermore, in the presence of high concentration of free 
carriers, as considered in this work, the excitonic effects 
get perfectly screened. The transition from bound ex- 
citonic state to free electron-hole plasma takes place at 
the Mott critical density^ If we apply its estimate from 
Ref. [H, we get about 5xl0 16 and 5xl0 18 cm -3 for InN 
and GaN, respectively. Above these values, excitons will 
dissociate which is easily met for the density range of this 
study. 

On the BGR, there are other theoretical estimations, 
such as the plasmon-pole^Z, and the Vashishta-Kalia 48 
models. However, we observed that when contrasted 
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with the existing experimental data on GaN 4 ^ and InNj 5 -^ 
the employed Berggren-Sernelius expression 34 in this and 
recent studies^ outperforms these other models in the 
relevant high density range of this study. From a for- 
mal point of view, there exist more rigorous many-body 
approaches which incorporate excitons, BGR and band- 
filling effects in the same level, however, their formidable 
computational cost leads to a trade off with a simplistic 
effective mass band structure^ 1 - 

Furthermore, we should mention that the employed 
Berggren-Sernelius expression assumes a cubic crystal 
with parabolic bands^ These simplifications have been 
avoided by Persson and co-workers who have incorpo- 
rated the electron-phonon contribution as well4^ Their 
framework is particularly better suited for handling the 
BGR under photoexcitation where both electron- and 
hole-gas many-body contributions can be individually in- 
cluded. However, as we stated previously, in the photoex- 
cited case we still use the Berggren-Sernelius expression 
but with the reduced electron- hole effective mass. Es- 
sentially, the same discussion also applies for the plasma 
contribution in Eq. (|5|). 

Compared to the available ab initio results for 
InN^ i 24 i 25 i 32 our band structure is in good overall agree- 
ment. Nevertheless, there remain further discrepancies 
of most of these theoretical results from the experimental 
data. First of all, due to the fact that our computational 
band gap of 0.85 eV for InN is somewhat higher than the 
experimentally established valued we expect especially 
our 1.55 fim wavelength data to be affected. However, 
this change should be more or less similar to the devia- 
tion between 1.55 /im and 1.3 /Ltm cases in Fig. 1131 In 
the case of GaN, our observation of a substantial devia- 
tion from the measured absorption in the vicinity of band 
gap (see, Fig. [3]) calls for improvement on pseudopoten- 
tial form factors for GaN which were originally optimized 
for high-field transport applications.— Furthermore, we 
do not include the band tail (Urbach's tail) absorption 
from the impurity states forming bands below the band 
edge . 52 i 53 Just like the excitons, the effect of localized 
states become marginal for high free-carrier densities, like 
beyond 10 18 cm~ 3 . Moreover, the characteristic Urbach 
energy parameter in high-quality InN samples turn out 
to be as low as 7.5 to 15 meV which indicates that these 
features are anchored to just below the band gap.— 

Being intended as a bulk study, this work neglects any 
surface effects, the most important of which is the sur- 
face charge accumulation causing built-in electric field 
and band bending . 55 ' 56 In this regard, it has been men- 
tioned by Kudrawiec et al. that for n-type InN layer 
since the region of strong band bending is the surface 
electron accumulation layer, it will not affect absorption- 
like techniques such as photoreflectance spectroscopy due 
to band filling of the conduction bandi 5 -^ Another effect 
discarded in this study is the thermal smearing of the de- 
generate carrier statistics. As a possible follow up of this 
work it can be readily included, nevertheless we do not 
expect any marked deviations. Even though the quasi- 



Fermi level shift in the valence band is of the same order 
as the thermal energy (cf., inset of Fig. [5]), as shown 
above, the main effect is driven by the Fermi level shift 
in the conduction band, which by far outweighs any ther- 
mal energy well up to room temperature for the carrier 
densities important for this work. Stemming from these 
idealizations, our main results in Fig. (TS] should be taken 
as upper bounds. The list of these simplifications also 
suggests possible improvements of this work. 

To harness the carrier-induced refractive index change 
such as in phase modulators, the associated losses need to 
be minimized. In this respect, the narrow 0.64 eV band 
gap of InN becomes undesirable for the 1.55 and 1.3 /jm 
applications as they enable supra band gap excitations 
where additional losses set in, especially for the photoex- 
cited case. Our comparative study here suggests that 
the In-rGai-zN can enjoy both the high carrier-induced 
refractive index sensitivity while at the same time operat- 
ing in the transparency region of the alloy. The poor sen- 
sitivity of GaN to carrier-induced refractive index change, 
as demonstrated by this work points towards the choice 
of indium-rich alloys. 



V. CONCLUSIONS 

We have shown that based on a full band theoreti- 
cal analysis, under n-type doping, InN and GaN should 
display different dielectric characteristics. The absorp- 
tion edge shift is masked in GaN, whereas it is highly 
pronounced in InN. For energies below 1 eV, the corre- 
sponding shifts in the real parts of the dielectric function 
for InN and GaN are in opposite directions. The free- 
carrier plasma contribution to refractive index becomes 
dominant above 10 20 cm -3 for the case of InN. At a wave- 
length of 1.55 fim, we predict more than 4% of change in 
the refractive index for InN for a doping of 10 19 cm~ 3 . 
Under optical pumping which fills the conduction and va- 
lence band edges with electrons and holes of equal den- 
sity, the refractive index is not further affected by the 
presence of the holes. However, its major consequence is 
the increased intravalence band aborption as a result of 
— > Tg transition. Beyond a hole density of 10 19 cm~ 3 , 
this affects 1.3 to 1.55 /Ltm wavelengths. We believe that 
this wide tunability of the index of refraction in InN by 
the conduction band electrons is technologically impor- 
tant for applications such as optical phase modulators. In 
the presence of holes, the valence band absorption needs 
to be taken into account, such as in the design of InN- 
based lasers for fiber-optics applications. Our findings 
suggest that the alloy composition of In^Gai-^N can be 
optimized in the indium-rich region so as to benefit from 
high carrier-induced refractive index change while oper- 
ating in the transparency region to minimize the losses. 
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APPENDIX: EMPIRICAL PSEUDOPOTENTIAL 
PARAMETERS FOR WURTZITE InN 

In the light of recent experimental and first-principles 
studies on wurtzite InN, we have expressed the local em- 
pirical pseudopotential method (EPM) form factors using 
the following functional form, 



TABLE I. Fitted values for the constants in the pseudopo- 
tential form factors V s and V a ; see Eqs. (1) and (2). Other 
parameters used in the EPM are also listed, where a and c are 
the lattice constants in the hexagonal plane and along the c 
axis, respectively, u is the wurtzite internal structural param- 
eter, and -Bmax denotes the radius of the energy sphere used 
for the reciprocal lattice vectors. 



Si 


= 0.261905 


ai 


= -0.411320 


S2 


= 0.232300 


ai 


= 0.176416 


S3 


= -2.726784 


az 


= 1.441420 


Si 


= 0.973335 


(14 


= 0.915637 


SB 


= 2.414422 


a B 


= 2.595997 


a - 


= 3.544 A 






c = 


= 5.716472 A 






u ■ 


= 0.379 






£ max = 12.5 Ry 







V s (q) = ( Sl q 3 + s 2 q 2 + S 3 ) exp(-s 4 q S5 ) , (A.7) 
V a (q) = (a iq 2 + a 2 q + 3 ) exp (-a 4 <f 5 ) , (A.8) 

where V s and V a are the symmetric and antisymmetric 
form factors in Rydbergs, q is the wave number in units 
of 2ir /a with a being the lattice constant in the hexagonal 
plane, and Si, a; (i = 1, . . . , 5) are the fitting parameters 
which are listed in Table I. 

The nonlocal part of the pseudopotential is based on 
the classical work of Chelikowsky and Cohen^ with the 
plane-wave matrix elements being given by 

V NL {K,K>) = —y2M(E)(2l + l)Pi(cos(9 KK ')) 

l.i 

xS\K - K')Ff(K, K'), (A.9) 



TABLE II. Fitted nonlocal EPM parameters for InN. A values 
are in Rydbergs and R values are in Angstroms. 



A^ n =5.55440, Rp n =0.195168 
A)f =4.51525, i^ n =0.132069 



A^=3.64897, R" =0.089251 



N_ 



where S l (q) is the structure factor™ for the atomic 
species i, and Pi{x) is the Legendre polynomial cor- 
responding to angular-momentum channel I, fl a is the 
atomic volume (i.e., primitive cell volume divided by 
number of atoms), K = k + G , K 1 = k + G , and 



Fi{K,K') 



[ii{KR)] 2 — ji-x(KR)ji + i(KR)^ 



K = K'. 



\Kj l+1 (KR) 3l {K'R) - K'ji +1 (K'R)ji (KR)} , K ± K< 



(A.10) 



here, ji(x) is the spherical Bessel function, R is (angular-momentum-dependent) well width. For InN, we treat the 
angular-momentum channels p and d for In, and p for N as nonlocal. Table II lists InN nonlocal EPM parameters. 
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